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We present a number of combinatorial characterizations of K-matrices. 
This extends a theorem of Fiedler and Ptak on linear-algebraic character- 
izations of K-matrices to the setting of oriented matroids. Our proof is 
elementary and simplifies the original proof substantially by exploiting the 
duality of oriented matroids. As an application, we show that a simple 
principal pivot method applied to the linear complementarity problems with 
K-matrices converges very quickly, by a purely combinatorial argument. 
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1 Introduction 

A matrix M € IR rixri is a P-matrix if all its principal minors (determinants of principal 
submatrices) are positive; it is a Z-matrix if all its off-diagonal elements are non-positive; 
and it is a K-matrix if it is both a P-matrix and a Z-matrix. 

Z- and K-matrices often occur in a wide variety of areas such as input-output produc- 
tion and growth models in economics, finite difference methods for partial differential 
equations, Markov processes in probability and statistics, and linear complementarity 
problems in operations research [2]. 

In 1962, Fiedler and Ptak [9] listed thirteen equivalent conditions for a Z-matrix to 
be a K-matrix. Some of them concern the sign structure of vectors: 
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Theorem 1.1 (Fiedler-Ptak [9]). Let M be a Z-matrix. Then the following conditions 
are equivalent: 

(a) There exists a vector x > such that Mx > 0; 

(b) there exists a vector x > such that Mx > 0; 

(c) the inverse M~ l exists and M _1 > 0; 

(d) for each vector i/O there exists an index k such that x^yk > for y = Mx; 

(e) all principal minors of M are positive (that is, M is a P-matrix, and thus a 
K-matrix). 

Our interest in K-matrices originates in the linear complementarity problem (LCP), 
which is for a given matrix M € IR nxn and a given vector q 6 R n to find two vectors w 
and z in W 1 so that 

w — Mz = q, 

w,z>0, (1) 
w T z = 0. 

In general, the problem to decide whether a LCP has a solution is NP-complete [6"lll5|. 
If the matrix M is a P-matrix, however, a unique solution to the LCP always exists |24j . 
Nevertheless, no polynomial-time algorithm to find it is known, nor are hardness results 
for this intriguing class of LCPs. It is unlikely to be NP-hard, because that would imply 
that NP = co-NP |18] . Recognizing whether a matrix is a P-matrix is co-NP-complete [8] . 
For some recent results, see also [2D]. 

If the matrix M is a Z-matrix, a polynomial-time (pivoting) algorithm exists [5] (see 
also [231 sect. 8.1]) that finds the solution or concludes that no solution exists. Interest- 
ingly, LCPs over this simple class of matrices have many practical applications (pricing 
of American options, portfolio selection problems, resource allocation problems). 

A frequently considered class of algorithms to solve LCPs is the class of simple principal 
pivoting methods (see Section [6] or Sect. 4.2]). We speak about a class of algorithms 
because the concrete algorithm depends on a chosen pivot rule. It has recently been 
proved in [11] that a simple principal pivoting method with any pivot rule takes at most 
a number of pivot steps linear in n to solve a LCP with a K-matrix M. 

The study of pivoting methods and pivot rules has led to the devising of combinatorial 
abstractions of LCPs. One such abstraction is unique-sink orientations of cubes [25]; the 
one we are concerned with here is one of oriented matroids. 

Oriented matroids were pioneered by Bland and Las Vergnas [I] and Folkman and 
Lawrence [10] , Todd [26] and Morris [19] gave a combinatorial generalization of LCPs 
by formulating the complementarity problem of oriented matroids (OMCP). Morris and 
Todd |21| studied properties of matroids extending LCPs with symmetric and positive 
definite matrices. Todd [26] proposed a generalization of Lemke's method [16] to solve 
the OMCP. Later Klafszky and Terlaky |14j and Fukuda and Terlaky [12] proposed a 
generalized criss-cross method; in [12] it is used for a constructive proof of a duality 
theorem for OMCPs in sufficient matroids (and hence also for LCPs with sufficient 
matrices). Hereby we revive their approach. 

In this paper, we present a combinatorial generalization (Theorem I5.4|) of the Fiedler- 
Ptak Theorem 11.11 To this end, we devise oriented-matroid counterparts of the con- 
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ditions (a)-(d). If the oriented matroid in question is realizable as the sign pattern of 
the null space of a matrix, then our conditions are equivalent to the conditions on the 
realizing matrix. In general, however, our theorem is stronger because it applies also to 
nonrealizable oriented matroids. 

As a by-product, our proof yields a new, purely combinatorial proof of Theorem I l.li 
Rather than on algebraic properties, it relies heavily on oriented matroid duality. 

We then use our characterization theorem to show that an OMCP on an n-dimensional 
K-matroid (that is, a matroid satisfying the equivalent conditions of Theorem 15. 4|) is 
solved by any pivoting method in at most 2n pivot steps. This implies the result of 
mentioned above that any simple principal pivoting method is fast for K-matrix LCPs. 

2 Oriented matroids 

The theory of oriented matroids provides a natural concept which not only generalizes 
combinatorial properties of many geometric configurations but presents itself in many 
other areas as well, such as topology and theoretical chemistry. 

2.1 Definitions and basic properties 

Here we state the definitions and basic properties of oriented matroids that we need in 
our exposition. For more on oriented matroids consult, for instance, the book [3]. 

Let £ be a finite set of size n. A sign vector on E is a vector X in {+1,0, —1} E . 
Instead of +1, we write just +; instead of —1, we write just — . We define X~ = 
{e G E : X e = -}, X G = {e G E : X e = - or X e = 0}, and the sets X°, X® and X + 
analogously. For any subset F of E we write Xp > if F C X® , and Xp < if 
F C X e ; furthermore X > if Xp > and X < if Xp < 0. The support of a sign 
vector X is X_ = X + UX~ . The opposite of X is the sign vector —X with (— X) + = X~ , 
(— X)~ = X + and (— X)° = X°. The composition of two sign vectors X and Y is given 



Definition 2.1. An oriented matroid on E is a pair M = (E,V), where V is a set of 
sign vectors on E satisfying the following axioms: 



by 




The product X ■ Y of two sign vectors is the sign vector given by 



{X ■ Y) e 



= < 



'0 if X e = or Y e = 0, 

+ if X e = Y e and X e / 0, 

— otherwise. 

V 



(VI) 0GV. 

(V2) If X G V, then -X G V. 
(V3) If X, Y G V, then X o Y G V. 
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(V4) If X,Y G V and e G X+ n Y~ , then there exists Z G V with Z+ C X+ U F+, 

z~ c rur, z e = o, and (x \y) u (y\x) u (i+nF + ) u (ruy - ) c z. 

The axioms (MU) up to ("V0J) are called vector axioms; (V[4|) is the vector elimination 
axiom. We say that the sign vector Z is the result of a vector elimination of X and 1" 
at element e. 

An important example is a matroid whose vectors are the sign vectors of elements of 
a vector subspace of W 1 . If A is an r x n real matrix, define 

V = {sgn x : x eR n and Ax = 0}, (2) 

where sgnx = (sgnxi, . . . , sgnx n ). Then V is the vector set of an oriented matroid on 
the set E = {1, 2, . . . , n}. In this case we speak of realizable oriented matroids. 

A circuit of A4 is a nonzero vector C G V such that there is no nonzero vector X G V 
satisfying IcC. 

Proposition 2.2. Lei = (E, V) 6e a matroid and let C be the collection of all its 
circuits. Then: 
(CI) 0C. 

(C2) IfCeC, then-C eC. 

(C3) For a// C,D eC, if C C L>, tten C = D or C = -D. 

(C4) If C,D e C, C ^ -D and e G C+ n £/ien ttere is a Z a C with Z + C 

(C+ U D+) \ {e} and Z~ C (C~ U D~) \ {e}. 
(C5) IfC,DeC,eeC + nD- and f G (C + \ D~) U (C~ \ i/ien i/tere is a Z a C 

with Z+ C (C+ U L>+) \ {e} ; Z- C(C"UD-)\ {e} and Z/ ^ 0. 
(C6) For every vector X G V f/iere exisi circuits C 1 ^ 2 , . . . ,C k G C snc/i i/iai X = 

C 1 o C 2 o • • • o C k and C\ • Cl > /or a// indices i,j and all e G X_. 

Moreover, if a set C of sign vectors on E satisfies ( (Uty-( (^), then it is the set of all 
circuits of a unique matroid; this matroid 's vectors are then all finite compositions of 
circuits from C. 

The property (C2J) is called weak circuit elimination; (C[5]) is called strong circuit 
elimination. In (C[6]) we speak about a conformal decomposition of a vector into circuits. 

A basis of an oriented matroid Ad is an inclusion-maximal set B C E for which there 
is no circuit C with CC B. Every basis B has the same size, called the rank of M. 

Proposition 2.3. Let B be a basis of an oriented matroid A4. For every e in E\B 
there is a unique circuit C(B,e) such that C(B,e) C B U {e} and C(B,e) e = +. 

Such a circuit C(B,e) is called the fundamental circuit of e with respect to B. 

Two sign vectors X and Y are orthogonal if the set {X e -Y e :e(z E} either equals {0} 
or contains both + and — . We then write X _L Y. 

Proposition 2.4. For every oriented matroid A4 = (E, V) of rank n there is a unique 
oriented matroid A4* = (E,V*) of rank \E\ — n given by 

V* = \y C {-,0, +} E : X 17 for every X G V j . 
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This M* is called the dual of M. Note that (M*)* = M. The circuits of M* are 
called the cocircuits of A4 and the vectors of At* are called the covectors of .M. The 
covectors of a realizable matroid given by ([2]) are sign vectors of the elements of the row 
space of the matrix A. 

We conclude this short overview by introducing the concept of matroid minors and 
extensions. For any FOE, the vector X \ F denotes the subvector (A e : e G E \ F) 
of X. Then let 

V \ F = {X \ F : X G V and X f = for all / G F} 
be the deletion and 

V/F = {X\F : X G V} 

the contraction of the vectors in V by the elements of F. It is easy to check that the 
pairs M. \ F = (E \ F,V \ F) and M/F = (E \ F,V / F) are oriented matroids. For any 
disjoint F,GCBwe call the oriented matroid {A4 \F)/G a minor of Ad. 

Note that M \ {e, e'} = (M \ {e}) \ {e'}, M/{e, e'} = (M/{e})/{e'} and (M \ {e})/ 
{e'} = (A4/{e'}) \ {e}, and so deletions and contractions can be performed element by 
element in any order, with the same result. 

Definition 2.5. A matroid A4 = (EU {q} , V) with q E is a one-point extension of KA 
ii M\{q} = M and there is a vector X of .M with X q ^ 0. 

2.2 Complementarity in oriented matroids 

In the rest of the paper, we are considering oriented matroids endowed with a special 
structure. The set of elements E<i n is a 2re-element set with a fixed partition E% n = SUT 
into two re-element sets and a mapping e4e from i?2n to Em which is an involution 
(that is, e = e for every e G Em) and for every e G 5 we have e G T. Note that this 
mapping constitutes a bijection between S and T. 

The element e is called the complement of e. For a subset F of i?2n let F = {e : e G -F}. 
A subset -F of E^n is called complementary if i* 1 n i 7 = 0. 

The matroids we are working with are of the kind A4 = (E^njV), where the set 
S C is a basis of A^. In addition, we study their one-point extensions M. = (E^n, V), 
where E2 n = -E^n U {q} for some element q ^ £?2n- Sometimes we make the canonical 
choice E2n = {!)•••) 2n} with 5" = {1, . . . , n} where the complement of an i G S is the 
element i + re. 

Definition 2.6. The oriented matroid complementarity problem ( OMCP) is to find a 
vector X of an oriented matroid A4 so that 

IeV, (3a) 
X > 0, X 9 = +, (3b) 
A e ■ Ae = for every e G -E^n, (3c) 

or to report that no such vector exists. 
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A vector X which satisfies ()3b|) is called feasible, one which satisfies ([3c]) is called 
complementary. Note that a vector is complementary if and only if its support is a 
complementary set. If an X € V satisfies (|3b|) and ([3c]) . then X is a solution to the 
OMCP(A4). 

Now we show how LCPs are special cases of OMCPs. Finding a solution to the 
LCP ([1]) is equivalent to finding an element x of 

V ={x£ R 2n+1 : [I n -M -q] x = o} 

such that 

x > 0, x 2n ,+i = 1, ,^ 
Xj • Xi +n = for every i € [n] . 

We set V = {sgnx : x € and consider the OMCP for the matroid Ai = (E2 n ,V). 
Clearly, if the OMCP has no solution, then V contains no vector x satisfying ([J]). If on 
the other hand there is a solution X satisfying ([3ap ~ ([3c]) . then the solution to the LCP 
can be obtained by solving the system of linear equations 

[I n -M -q]x = 0, 

Xi = whenever = 0, 

X2n+l = 1- 

This correspondence motivates the following definition. 

Definition 2.7. An oriented matroid Ai = (B2„,V) is L CP-realizable if there is a matrix 
Mel m such that 

V = {sgnx : x £ M 2n and [l n -Af] x = oj. 

The matrix M is then a realization matrix of M.. This is a little nonstandard, because 
usually the matrix A from ([2]) is called a realization matrix. The columns of I n are 
indexed by the elements of S C E%ni an d the columns of — M are indexed by the elements 
of T C so that if the fcth column of I ra is indexed by e, then the kth column of — M 
is indexed by e. 

The extension = (i?2n)V) is L CP-realizable if there is a matrix M € M nxn and a 
vector geM" such that 

V = {sgnx : x € M 2n+1 and [/„ -M -q]x = o\. 

To study the algorithmic complexity of OMCPs, we must specify how the matroid M. 
is made available to the algorithm. We will assume that it is given by an oracle which, for 
a basis B of Ai and a nonbasic element e G Z?2n \ B, outputs the unique (fundamental) 
circuit C of M. with support CCBU {e} such that X e = +. 

In the LCP-realizable case such an oracle can be implemented in polynomial time; 
in fact, it consists in solving a system of 0{n) linear equations in 2n + 1 variables. 



6 



Thus, in the RAM model, the oracle can be implemented so that it performs arithmetic 
operations whose number is bounded by a polynomial in n. Hence our goal is to develop 
an algorithm that solves an OMCP using a number of queries to the oracle that is 
polynomial in n. 

Such an algorithm for the OMCP would obviously provide a strongly polynomial 
algorithm for the LCP. Since the LCP is NP-hard in general, the existence of such an 
algorithm is unlikely. In Section [6] we do, nevertheless, prove the existence of such an 
algorithm for a special class of oriented matroids: K-matroids. 

3 P-matroids 

In this and the following sections, we investigate what properties of oriented matroids 
characterize oriented matroids realizable by special classes of matrices. We start with 
P-matrices; recall that a P-matrix is a matrix whose principal minors are all positive. 
Several conditions are equivalent to the positivity of principal minors: 

Theorem 3.1. For a matrix M £ M nxn , the following are equivalent: 

(a) All principal minors of M are positive (i.e., M is a P-matrix); 

(b) there is no nonzero vector x such that x^Uk < for all i = 1,2, ... ,n, where 
y = Mx; 

(c) the LCP (H|) with the matrix M and any right-hand side q has exactly one solution. 

The equivalence of (jaj) and (jb]) is due to Fiedler and Ptak |9j. The equivalence of (jaj) 
and (jcj) was proved independently by Samelson, Thrall and Wesler [21], Ingleton [T5] , 
and Murty [22]. 

The following notions and our definition of a P-matroid are motivated by the condi- 
tion (jb]) in Theorem 13.11 It is much easier to express in the oriented-matroid language 
than (jaj). 

A sign vector X G { — ,0, +} E2n is sign-reversing (s.r.) if X e ■ Xe < for every e € S. 
If in addition X_ = Em, the vector is totally sign-reversing (t.s.r.). Analogously, an X 
is sign-preserving (s.p.) if X e ■ Xe > for every e, and totally sign-preserving (t.s.p.) if 
X = Em as well. 

Definition 3.2 (Todd |26j). An oriented matroid M. on -E^n is a P-matroid if it has no 
sign-reversing circuit. 

Note that a P-matroid contains no nonzero sign-reversing vectors, because every vector 
is the composition of some circuits and composing non-s.r. circuits yields non-s.r. vectors. 
Hence, using Theorem 13. lj we immediately get: 

Proposition 3.3. 

(i) If M is L CP-realizable and there exists a realization matrix M that is a P-matrix, 
then M is a P-matroid. 

(ii) If A4 is an L CP-realizable P-matroid, then every realization matrix M is a P- 
matrix. 
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P-matroids were extensively studied by Todd [26]. He lists eight equivalent conditions 
for a matroid to be a P-matroid. We recall three of them (conditions (a), (a*) and (c) 
below) and add two new ones. 

Theorem 3.4. For an oriented matroid M. on Em, the following conditions are equiv- 
alent: 

(a) M. has no s.r. circuit; 
(a*) M. has no s.p. cocircuit; 

(b) every t.s.p. X is a vector of KA; 
(b*) every t.s.r. Y is a covector of KA; 

(c) every one-point extension A4 of M to E^n contains exactly one complementary 
circuit C such that C > and C q = +. 

Proof. The equivalence of the conditions (a), (a*) and (c) was shown by Todd [26j . 
Morris [19] proved that (a) implies (b). We show the equivalence of (a) with (b*). The 
equivalence of (a*) with (b) is proved analogously. 

First we prove that (a) implies (b*). Since no circuit of Ad is s.r., there is for every 
circuit C an element e such that C e ■ CV = +. It follows that every t.s.r. sign vector Y 
is orthogonal to every circuit, hence Y is a covector. 

For the opposite direction, suppose that there is a s.r. circuit C. If so, then any 
t.s.r. vector Y for which Y + C C + and Y~ C C~ is not orthogonal to C, which is a 
contradiction with (b*). □ 

The condition (b) of this theorem has a translation for realization matrices of P-ma- 
troids, that is, for P-matrices: 

Corollary 3.5. A matrix M € M. nxn is a P-matrix if and only if for every a £ { — 1, +l} n 
there exists a vector x € R™ such that for y = Mx and for each i E {1, 2, . . . , n} we have 

uiXi > 0, 
nVi > 0. 

Todd [26] also gives an oriented-matroid analogue of the "positive principal minors" 
condition. Stating it would require some more explanations; later in this article we need 
a weaker property of P-matroids, though, which corresponds to the fact that all principal 
minors of a P-matrix are nonzero. 

Lemma 3.6 (Todd [26]). For a P-matroid M. every complementary subset B C of 
cardinality n is a basis. 

Remark. In addition, every such complementary B is also a cobasis, i.e., it is a basis 
of the dual matroid M* . 

Next we consider principal pivot transforms (see [271 [28]) °f P-matrices. The fact that 
every principal pivot transform of a P-matrix is again a P-matrix [2U] is well-known. 
The proof is not very difficult but it uses involved properties of the Schur complement. 
In the setting of oriented matroids the equivalent is much simpler. First let us define 
principal pivot transforms of oriented matroids. 
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Definition 3.7. Let F C Em be a complementary set. The principal pivot transform 
of a sign vector X with respect to F is the sign vector C given by 



The principal pivot transform of a matroid ftA with respect to F is the matroid whose 
circuits (vectors) are the principal pivot transforms of the circuits (vectors) of A4. 

It is easy to see that, in the LCP-realizable case, principal pivot transforms of a 
matroid correspond to matroids realized by corresponding principal pivot transforms of 
the realization matrix. Thus the following proposition implies the analogous theorem 
for P-matrices. 

Proposition 3.8. Every principal pivot transform of a P-matroid is a P-matroid. 

Proof. The principal pivot transform of a circuit C is sign-reversing if and only if C is 
sign-reversing. □ 

4 Z-matroids 

The second class of matrices we examine are Z-matrices; the corresponding matroid 
generalizations are Z-matroids. Recall that a Z-matrix is a matrix whose every off- 
diagonal element is non-positive. The definition of Z-matroids was first proposed in |17| . 

Definition 4.1. A matroid A4 on E2 n is a Z-matroid if for every circuit C of KA we 
have: 



Remark. In the definition of Z-matroid we could replace all occurrences of the word 
"circuit" with the word "vector". Indeed, in a conformal decomposition of a vector 
violating ([5]), there would always be a circuit violating ([5]) as well. 

It makes perfect sense to define Z-matroids in this way. We show that in LCP- 
realizable cases, any realization matrix M is a Z-matrix. 

Proposition 4.2. 

(i) If M. is LCP -realizable and there exists a realization matrix M that is a Z-matrix, 
then M is a Z-matroid. 

(ii) If Ad is an LCP-realizable Z-matroid, then every realization matrix M is a Z- 



Proof. We fix Em = {!)•■•) 2n} with S = {1, . . . , n} where the complement of an i G S 
is the element i + n. 




If Or > 0, then 



(5) 



Ce = + for all e G S with C e = +. 



matrix. 
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(i) Let ei denote the ith unit vector and m,j the jth column of the matrix M. The 
sign pattern of the Z-matrix M implies that there is no linear combination of the form 



2)i. 



+ J2 x i e i ~ Yl x o m 3 =0 ' 



3=1 j=n+l 

where Xj > for every j > n and Xj+ n = 0, because the ith row of the left-hand side 
is strictly positive. Hence there is no vector X £ V for which Xt > 0, = + but 
Xj +n = for some i £ S. 

(ii) Proof by contradiction. Assume that for an LCP-realizable Z-matroid M (where 
5 is a basis), there is a realization matrix M that is not a Z-matrix, that is, there is an 
off-diagonal rriij > 0. If so, there is a vector X with Xj +n = + and X T \ij +n \ = 0, but 
Xi = +. This X violates the Z-matroid property ([SJ since also Xi +n = 0, a contradiction. 
Thus no positive m,j can exist and M has to be a Z-matrix. □ 

Another option is to characterize a Z-matroid with respect to the dual matroid M*. 

Proposition 4.3. An oriented matroid M on E2 n is a Z-matroid if and only if for every 
cocircuit D of M we have: 

If D s < 0, then 

(6) 

Dp = — for all e € T with D e = +. 

Proof. First we prove the "only if" direction. Suppose that there is a cocircuit D which 
does not satisfy ©. Accordingly Ds < and there is e € T such that D e = +, but 
= 0. But note that the fundamental circuit C := C(S,e) and D are not orthogonal 
because the Z-matroid property ([SJ implies that Cs\{e} — 0. Hence no such D can exist. 

For the "if" direction suppose that there is a circuit C for which Ct > and C e = +, 
but CV = for some e £ S. This circuit C and the fundamental cocircuit D := D(T,e) 
are not orthogonal since by assumption (J5J) holds for D and of course —D, hence Df\{e} > 
0. □ 

In the proofs in the following section we often make use of fundamental circuits. Here 
we observe that all fundamental circuits with respect to the basis S follow the same sign 
pattern. 

Lemma 4.4. Let M. be a Z-matroid. Let e € T and let C = C(S,e) be the fundamental 
circuit of e with respect to the basis S. Then 

C e = +, 

Cr\{ e } = 0, 

C S \{e} < 0. 

Proof. The first and the second equality follow directly from the definition of a funda- 
mental circuit. Thus Ct > 0. Hence the third property follows from the Z-matroid 
property ©. □ 
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5 K-matroids 



Definition 5.1. A matroid A4 on E^n is a K-matroid if it is a P-matroid and a Z- 
matroid. 

Combining Proposition 13.31 and Proposition 14,21 we immediately get: 
Proposition 5.2. 

(i) If M is I CP-realizable and there exists a realization matrix M that is a K-matrix, 
then M is a K-matroid. 

(ii) If M is an L CP-realizable K-matroid, then every realization matrix M is a K- 
matrix. 

An oriented matroid minor Ai \ F/F where F is a complementary subset of E^n is 
called a principal minor of M. 

Lemma 5.3. Let A4 be a K-matroid. Then every principal minor A4 \ F/F is a K- 
matroid. 

Proof. It was shown by Todd [26] that every principal minor of a P-matroid is a P- 
matroid. Thus, it is enough to show that such a minor is a Z-matroid, and for this, since 
deletions and contractions can be carried out element by element in any order, it suffices 
to consider the case that F is a singleton. 

First, we prove that if e G T, then A4 \ {e} / {e} is a Z-matroid. Such a principal 
minor consists of all circuits C \ {e, e}, where C is a circuit of M and C e = 0. Since 
every circuit of M satisfies the Z-matroid characterization ([5]), such a circuit C \ {e, e} 
trivially satisfies it too. 

Secondly, let e G S. Here we apply a case distinction. If Ce = +, then (C\{e,e})r > 
if and only if Ct > 0. As a direct consequence, C \ {e,e} satisfies ([5]) because C does. 
If Ce = — , we can show that there is another element / G T such that Cj = — too, 
that is, (C \ {e,e})T ^ and thus the Z-matroid property ([5]) is obviously satisfied. 
Assume for the sake of contradiction that there is no such / £ T. The strong circuit 
elimination (C[5]) of C and the fundamental circuit C(S,e) at e then yields a circuit C 
with C' T > 0, C^ = and Cg = +. Since e £ S, such a C would violate the Z-matroid 
definition, a contradiction. □ 

Our main result, the combinatorial generalization of the Fiedler-Ptak Theorem II .11 is 
the following. 

Theorem 5.4. For a Z-matroid A4 (with vectors V, covectors V* , circuits C and cocir- 
cuits T> ), the following statements are equivalent: 

(a) 3X £ V : Xt > and X s > 0; (a*) 3Y € V* : Y s < and Y T > 0; 

(b) 3X G V : X > 0; (b*) 3Y G V* : Y s < and Y T > 0; 

(c) VC G C: Cs > =^ C T > 0; (c*) V£> G V: D T > =^ £> s < 0; 

(d) there is no s.r. circuit C G C (d*) i/iere is no s.p. cocircuit D G T>. 
(that is, M. is a P-matroid); 
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In order to use duality in the proof of this theorem, let us first define the reflection of 
a matroid M = (E 2n , V) to be the matroid &(M) = {E 2n , 3?(V)), where K(V) = : 
IeV} with 

v v Ve \-Xe ifeer. 

Observe that = because of (Ml, and that *d(M*) = ^(M)*; thus 

$($l(M*)*) = M. (7) 

Proof of Theorem \5.4\ 

(a) =>■ (b): Let X be as in (a). Since Xt > 0, the Z-matroid property ([5]) implies 
that if X e = + for aneeS, then X F = +. Thus X T > 0. 

(b) => (c): Let X be the all-plus vector as in (b). Suppose that there is a circuit 
C G C not satisfying (c), that is, Cs > but C e = — for some element e in T. 
Starting with 1"° = C, we apply a sequence of vector eliminations (M3) to get 
vectors Y l . We eliminate Y l_1 and X at any element e where Yl~ x = — . For a 
resulting vector Y J it holds that (Y l )~ C (Y 1 ^ 1 )^ . Thus, at some point (Y k )~ = 
while Y k = and Y^ = + where e € T is the element eliminated in step k — 1. This 
vector y fc does not satisfy the Z-matroid property ([5]), which is a contradiction. 

(c) =>■ (d): Suppose that there is a s.r. circuit C G C, that is, C e • Ce < for every 
e € 5. Let C° = C. We apply consecutive circuit eliminations (CHJ). To get C\ we 
eliminate C l_1 with any fundamental circuit C := C(S,e) at position e € T where 
Cg _1 = -. By Lemma PI we have Cg\{e} < 0. 

After finitely many eliminations we end up with a circuit C k for which Cj> > 0. 
Now we claim that Cg < 0: Indeed, if e €. S such that = +, then Ce = +, 
and thus C e < because C is sign-reversing. Since we never eliminate at e, all 
fundamental circuits C used in the eliminations satisfy C e < as noted above. 
Hence C k e < 0. If on the other hand C| = 0, then C e fc < by ©. 

Moreover, since S is a basis, C^_ ^ 5, and so there exists e G T with Cg = +. 
Therefore — C k violates property (c), a contradiction. 

(d) =>■ (a*): Because of (d), for every circuit C there is an e G S such that C e -Ce = +• 
The sign vector K where Y$ < and It > is orthogonal to every circuit C, 
because the sign of Y e ■ C e is opposite to the sign of Yg ■ Ce- Hence such a Y is a 
covector. 



To finish the proof, notice that a matroid A4 satisfies (a*) if and only if the reflection 
of its dual 3t(M*) satisfies (a); analogously for (b*) and (b), (c*) and (c), and (d*) 
and (d). Thus if M. satisfies (a*), then $t(A4*) satisfies (a), hence also (b), and so 
(using 0) M satisfies (b*). The missing implications (b*) => (c*), (c*) (d*), 
and (d*) =^ (a) are proved analogously. □ 
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6 Algorithmic aspects 



Let an OMCP(TW) be given, where A4 is any one-element extension of an n-dimensional 
matroid M. on Em- We present simple principal pivot algorithms to find a solution. 
This kind of algorithm is a well-established solving method for LCPs. Sometimes called 
Bard-type methods, they were first studied by Zoutendijk [30J and Bard [lj. 

Here we extend a recent result of [TT] to the generalizing setting of OMCP. We prove 
below that the unique solution to every OMCP(A-'f) where the underlying matroid A4 
is a K-matroid, is found by every simple principal pivot algorithm in a linear number of 
pivot steps. 

Let Ai be given by an oracle which, for a basis B of A4 and a non-basic element 
e £ E2n\B, outputs the unique fundamental circuit C(B,e). A simple principal pivot al- 
gorithm starts with a fundamental circuit C° = C(B°, q) where B° is any complementary 
basis. For instance, we can take B = S. It then proceeds in pivot steps. Assume that 
the ith. step leads to a fundamental circuit C l = C(B l , q). We require the complementary 
condition ([3c]) to be an invariant, that is, B % is supposed to be complementary. If C l is 
feasible, that is, the condition (j3bj) is satisfied, then C % is the solution and the algorithm 
terminates. Otherwise, we obtain C t+l as follows: choose an e* £ B % for which C* ; = — 
according to a pivot rule. Then the pivot element e* is replaced in the basis with its 
complement e l , that is, B l+1 = B % \ {e*} U {e 4 }. Lemma 13.61 asserts that B l+1 is indeed 
a basis. Then C l+l = C(B 1+1 ,q) is computed by feeding the oracle with basis B t+l and 
the non-basic element q. The algorithm then proceeds with pivot step i + 2. 

SimplePrincipalPivot(A4, B°) 
i := 

C° := C{B°,q) 
while {C i y / do 

choose e l E (C l )~ according to a pivot rule R 

B i+i := B i \ (e^ U {e?} 

C i+1 := C(B i+ \q) 

i:=i + l 
end while 
return C % 

If the number of pivots is polynomial in n, then the whole algorithm runs in polynomial 
time too, provided that the oracle computes the fundamental circuit in polynomial time. 
This is the case if the LCP is given by a matrix M and a right-hand side q as in ([1]). 

The number of pivots depends on the applied pivot rule and some rules may even 
enter a loop on some inputs A4. If the input is a K-matroid extension, though, then 
the SimplePrincipalPivot method is fast. We claim that no matter which pivot rule 
is applied, SimplePrincipalPivot runs in a linear number of pivot steps on every K- 
matroid extension. The following two lemmas are required to prove this fact. While 
the first one holds for every P-matroid extension, the second is restricted to K-matroid 
extensions. 
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Lemma 6.1. If M is a P-matroid extension, then C^t 1 = + for every i > 0. 

Proof. First suppose that C^t 1 = — in some pivot step i + 1. Let C" be the result of a 

weak circuit elimination of C l and —C l+1 at q. Then C is contained in B l U {e*}, and 
C'i < and C— > 0, in other words it is a s.r. vector. According to the Definition 13,21 

of a P-matroid, no s.r. circuit can exist. Thus > 0. 

Now suppose that CLf 1 = 0. Then C %+1 is also the fundamental circuit C(B l ,q), 

hence C t+1 = C\ This is a contradiction because C 4 ^ 1 = ^ - = C\. □ 

Lemma 6.2. // A4 is a K-matroid extension, then for every f £ T: 

IfCt>0 for some h>0, then Cj > /or every > /j. 

Proof. For the sake of contradiction suppose that the statement does not hold and let 
I > h be the smallest value such that Cj- > 0, but Cy +1 = — . By Lemma IBTTl f ^ e l , and 
so / lies in B l and in Let X be the result of a vector elimination of C and — C' +1 

at g. Note that X e i = — , Xf = + and Xj = 0. In addition by Lemma 16. II it holds that 
X-j = —. Since X q = 0, the sign vector X \ {q} is a vector of the K-matroid AL Now 

let F := B l \ {/, e 1 }. Consider the principal minor Ai \ F/F, which is a matroid on the 
element set {/, /, e , e'}. By Lemma 15.31 it is also a K-matroid. Further it contains the 
vector X' = X \ [e 2ti \ {/,/,e*,?}) with X' &1 = -, X' f = +, XL = and X<- = -. 

The contradiction follows from the fact that —X' violates the K-matroid property (c) in 
Theorem 15.41 □ 

Theorem 6.3. Every simple principal pivot algorithm runs in at most In pivot steps 
on every K-matroid extension. 

Proof. We prove that, no matter which pivot rule R one applies, every element e G 
E2n is chosen at most once as the pivot element. Consider any pivot step h in the 
SimplePrincipalPivot algorithm. First suppose that the pivot element e h is in S. 
According to Lemma 16.11 C^ 1 > 0. Moreover, for every k > h we have CK > 

(Lemma I6.2p and C^ h = 0. In other words, the elements e h and e h cannot become pivot 
elements in later steps. Second, suppose that the pivot e h is in T. Then the argument 
above fails. Even though C-^" 1 > (Lemma 16. lh . we cannot conclude that C^- > for 
every k > h, because Lemma 16.21 does not apply. It may eventually happen for some k 
that e h is chosen as pivot e k . However if so, our first argument applies for pivot step k 
and neither e h nor e h can become pivot elements again. □ 

Remark. If SimplePrincipalPivot starts with the basis B° = S, then at most n 
pivot steps are needed, because Cj. = and hence, by Lemma HO} C % T > for all i. 
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7 Extension to principal pivot closures 



So far, we have considered a matroid M on a complementary set E^n where the maximal 
complementary set S is fixed from the beginning. In the following, S' is an arbitrary 
complementary subset of size n and T' = {e : e G S"}. 

Definition 7.1. A matroid .M on i?2n is a Z* -matroid if there is a complementary set 
S' C i? 2n of cardinality n such that for T" = {e : e G S"} and every circuit C of Al we 
have: 

If C T > > 0, then 

CV = + for all e G S' with C e = +. 

Analogously A4 is a K*-matroid if it is a P-matroid and a Z*-matroid. Note that the 
class of Z*- and K*-matroids are the closures, under all principal pivot transforms, of Z- 
and K-matroids, respectively. Moreover, Proposition 14. 3| Lemma [5,3l up to Theorem 16.31 
have equivalent counterparts for these closure classes, obtained by substituting S by S' 
and accordingly T by T' in the original statements. Hence we get the following. 

Corollary 7.2. Every simple principal pivot algorithm finds the solution to OMCP(M. ), 
where M is a K*-matroid extension, in at most In pivot steps. 

The reader might wonder why we introduced Z-matroids and K-matroids at all and 
did not start off with their principal pivot closures. One good reason for our approach 
is to point out the correspondence of LCP-realizable Z-matroids and their matrix coun- 
terparts, see Proposition 14.21 With respect to this, the main problem is that a principal 
pivot transform of a Z-matroid or a K-matroid is in general not a Z-matroid, a K-matroid 
respectively. However every principal pivot algorithm is still able to solve an LCP(M, q) 
where M is a principal pivot transform of a K-matrix in a linear number of pivot steps. 

References 

[1] Yonathan Bard. An eclectic approach to nonlinear programming. In Optimization 
(Proc. Sem., Austral. Nat. Univ., Canberra, 1971), pages 116-128, St. Lucia, 1972. 
Univ. Queensland Press. 

[2] Abraham Berman and Robert J. Plemmons. Nonnegative matrices in the mathe- 
matical sciences. Academic Press, 1979. 

[3] Anders Bjorner, Michel Las Vergnas, Bernd Sturmfels, Neil White, and Giinter M. 
Ziegler. Oriented Matroids, volume 46 of Encyclopedia of mathematics and its ap- 
plications. Cambridge University Press, second edition, 1999. 

[4] Robert G. Bland and Michel Las Vergnas. Orientability of matroids. J. Combin. 
Theory, Ser. B, 24:94-123, 1978. 



15 



[5] Ramaswamy Chandrasekaran. A special case of the complementary pivot problem. 
Opsearch, 7:263-268, 1970. 

[6] Sung-Jin Chung. NP-completeness of the linear complementarity problem. J. Op- 
tim. Theory Appl, 60(3):393-399, 1989. 

[7] Richard W. Cottle, Jong-Shi Pang, and Richard E. Stone. The Linear Comple- 
mentarity Problem. Computer science and scientific computing. Academic Press, 
1992. 

[8] Gregory E. Coxson. The P-matrix problem is co-NP-complete. Math. Programming, 
Ser. A, 64(2): 173-178, 1994. 

[9] Miroslav Fiedler and Vlastimil Ptak. On matrices with non-positive off-diagonal 
elements and positive principal minors. Czechoslovak Math. J., 12 (87):382-400, 
1962. 

[10] Jon Folkman and Jim Lawrence. Oriented matroids. J. Combin. Theory, Ser. B, 
25:199-236, 1978. 

[11] Jan Foniok, Komei Fukuda, Bernd Gartner, and Hans- Jakob Liithi. Pivoting in 
linear complementarity: Two polynomial-time cases. Discrete Comput. Geom., 
42(2):187-205, 2009. 

[12] Komei Fukuda and Tamas Terlaky. Linear complementarity and oriented matroids. 
J. Oper. Res. Soc. Japan, 35(1):45-61, 1992. 

[13] Aubrey William Ingleton. A problem in linear inequalities. Proc. London Math. 
Soc. (3), 16:519-536, 1966. 

[14] Emil Klafszky and Tamas Terlaky. Some generalizations of the criss-cross method for 
the linear complementarity problem of oriented matroids. Combinatorica, 9(2):189- 
198, 1989. 

[15] Masakazu Kojima, Nimrod Megiddo, Toshihito Noma, and Akiko Yoshise. A unified 
approach to interior point algorithms for linear complementarity problems, volume 
538 of Lecture Notes in Computer Science. Springer- Verlag, Berlin, 1991. 

[16] Carlton E. Lemke. Bimatrix equilibrium points and mathematical programming. 
Management Sci., 11:681-689, 1964/1965. 

[17] Carlton E. Lemke and Hans-Jakob Liithi. Classes of matroids and complementarity. 
30 pages manuscript, February 1986. 

[18] Nimrod Megiddo. A note on the complexity of P-matrix LCP and computing an 
equilibrium. RJ 6439, IBM Research, Almaden Research Center, 650 Harry Road, 
San Jose, California, 1988. 



16 



[19] Walter D. Morris, Jr. Oriented Matroids and the Linear Complementarity Problem. 
PhD thesis, Cornell University, Ithaca, 1986. 

[20] Walter D. Morris, Jr. and Makoto Namiki. Good hidden P-matrix sandwiches. 
Linear Algebra AppL, 426(2-3) :325-341, 2007. 

[21] Walter D. Morris, Jr. and Michael J. Todd. Symmetry and positive definiteness in 
oriented matroids. European J. Combin., 9(2): 121-129, 1988. 

[22] Katta G. Murty. On the number of solutions to the complementarity problem and 
spanning properties of complementary cones. Linear Algebra and AppL, 5:65-108, 
1972. 

[23] Katta G. Murty. Linear Complementarity, Linear and Nonlinear Programming, 
volume 3 of Sigma Series in Applied Mathematics. Heldermann, Berlin, 1988. 

[24] Hans Samelson, Robert M. Thrall, and Oscar Wesler. A partition theorem for 
Euclidean n-space. Proc. Amer. Math. Soc, 9:805-807, 1958. 

[25] Alan Stickney and Layne Watson. Digraph models of Bard-type algorithms for the 
linear complementarity problem. Math. Oper. Res., 3(4):322-333, 1978. 

[26] Michael J. Todd. Complementarity in oriented matroids. SIAM J. Alg. Disc. Meth., 
5:467-485, 1984. 

[27] Michael J. Tsatsomeros. Principal pivot transforms: Properties and applications. 
Linear Algebra AppL, 307(1-3):151-165, 2000. 

[28] Albert W. Tucker. A combinatorial equivalence of matrices. In Richard Bellman and 
Marshall Hall, Jr., editors, Proc. Sympos. Appl. Math., volume 10, pages 129-140. 
American Mathematical Society, 1960. 

[29] Albert W. Tucker. Principal pivotal transforms of square matrices. 1963 Spring 
Meeting of SIAM, Menlo Park, California, 1963. Abstract in SIAM Review, vol. 5, 
p. 305. 

[30] Guus Zoutendijk. Methods of feasible directions: A study in linear and non-linear 
programming. Elsevier Publishing Co., Amsterdam, 1960. 



17 



